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We implement self-consistent microscopic calculations in order to describe out-of-equilibrium non- 
local transport in normal metal-superconductor-normal metal hybrid structures in the presence of 
a magnetic field and for arbitrary interface transparencies. A four terminal setup simulating usual 
experimental situations is described by means of a tight-binding model. We present results for 
the self-consistent order parameter and current profiles within the sample. These profiles illustrate 
a crossover from a quasi-equilibrium to a strong non-equilibrium situation when increasing the 
interface transparencies and the applied voltages. We analyze in detail the behavior of the non- 
local conductance in these two different regimes. While in quasi-equilibrium conditions this can be 
expressed as the difference between elastic cotunneling and crossed Andreev transmission coefficients, 
in a general situation additional contributions due to the voltage dependence of the self-consistent 
order parameter have to be taken into account. The present results provide a first step towards a self- 
consistent theory of non-local transport including non-equilibrium effects and describe qualitatively 
a recent experiment [Phys. Rev. Lett. 97, 237003 (2006)]. 



I. INTRODUCTION 

In addition to charge and spin, future electronic de- 
vices may also manipulate the non-local correlations al- 
lowed by quantum mechanics, known as entanglement. 
For instance a superconductor may be used as a source 
of Einstein Podolsky Rosen pairs of electronsi*^. Local 
Andreev reflection^ at a normal metal - superconductor 
(NS) interface is a process by which a spin-up electron 
incoming from the normal side is reflected as a spin-down 
hole while a pair is transferred into the superconductor: 
pairs of electrons penetrate the superconductor for an ap- 
plied bias smaller than the gap (see Fig.[T^). For an op- 
posite applied voltage, the superconductor S emits corre- 
lated pairs of electrons into the normal metal N. Andreev 
reflection takes place in a coherence volume of size i the 
characteristic length associated to the superconducting 
gap A. Two separate normal electrodes connected to a 
superconductor within a distance of order niay thus be 
coupled by "non-local" or "crossed" Andreev processes 
1,2,4-30 (CAR, see Fig. [lb). On the other hand another 
type of non-local process may take place. Electrons can 
also tunnel across the superconductor from one normal 
electrode to the other. This normal tunneling has been 
called^ "elastic cotunneling" (EC) in analogy to similar 
processes taking place in Coulomb blockaded quantum 
dots^. 

Motivated by the possibility of creation of non-local 
coherent states by means of CAR processes, several ex- 
periments have been performed recently^i^i^ on supercon- 
ducting structures connected to normal or ferromagnetic 



metallic electrodes. Depending on the type of S/N in- 
terfaces, geometry and range of parameters, different be- 
haviors of the non-local conductance or resistance have 
been observed. 

According to lowest order perturbation theory^ in the 
tunnel amplitudes, EC and CAR have the same trans- 
mission coefficient, once an average over the Fermi wave- 
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FIG. 1: (Color online). Schematic representation of the elec- 
trical circuit used in a non-local conductance experiment. The 
lowest order processes are shown on the figure: local Andreev 
refiection (a), crossed Andreev refiection (CAR) (b) and elas- 
tic cotunneling (EC) (c). 
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length scale or over disorder is carried out. Since an 
opposite charge is transmitted by EC and CAR, it is de- 
duced that the non-local conductance vanishes in this 
limit. A finite non-local signal may be restored by 
higher order tunneling processe a^i^i^^ or by Coulomb 
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interactionsii. Moreover, out-of-equilibrium effects may 
play an important role on the non-local transport as was 
suggested in Ref. [13 for a superconducting quantum dot 
and in Ref. [l^ for normal electrodes connected to a three 
dimensional superconductor. 

Experimental data on the other hand arc still not well 
understood. In particular the experiment of Ref. on a 
planar NISIN structure (where I stands for an insulat- 
ing barrier) has provided unexpected experimental ev- 
idence for a non-local signal dominated either by EC 
or by CAR depending on the value of the applied bias. 
Also importantly, this experimenij^ has shown a suppres- 
sion of the non-local signal when an external magnetic 
field, much smaller than the critical one, was applied. 
Other experiments on NSN hybrid structures provide ev- 
idence for charge imbalance effects as the voltage ap- 
proaches the temperature-dependent gap valued and in 
a temperature window close to the superconducting crit- 
ical temperature^. 

Although non-local transport in SN structures has 
been addressed theoretically by several works in recent 
years, none of these provide a full self-consistent model 
for describing the superconducting order parameter, the 
current profile and the non-local conductance of the S 
region, as well as the magnetic field dependence of these 
quantities. 

To bridge this gap. we analyze here a NISIN planar 
structure using a microscopic model which we solve self- 
consistently in order to obtain the current profile inside 
the superconductor in the presence of an applied volt- 
age and magnetic field. To adequately represent a non- 
local transport measurement setup in the self-consistent 
calculation it is important to include additional super- 
conducting leads which allow to remove from the sample 
the injected current (see Fig. [2]). This implies a substan- 
tial difference with conventional "two terminal" measure- 
ments, for which the current is the same on both leads. A 
two-terminal situation would be recovered by removing 
the superconducting leads. 

An important parameter characterizing an actual ex- 
perimental situation is the value of the transparency for 
the barriers connecting the superconductor with the nor- 
mal electrodes. Although this parameter is difficult to 
be controlled in experiments the available data corre- 
spond to very different ranges as it was pointed out in 
Ref. [sO- Within our model it is possible to study non- 
local transport for arbitrary interface transparency, ap- 
plied bias and external magnetic field. This allows us to 
analyze the crossover from a quasi-equilibrium situation 
at low barrier transparency to the case of strong non- 
equilibrium for high transparency and finite voltages. In 
agreement with previous works^i^, we find that increas- 
ing barrier transparency leads to a dominance of EC over 




FIG. 2: (Color online) . Schematic representation of the tight- 
binding model used in our simulations and of the electrical 
circuit (green lines). The self-consistent gap is evaluated on 
a square lattice with M — Ry/ao transverse channels and 
a length R^c ~ Nao, with ao the lattice spacing. Normal 
electrodes Nj, (Na) are connected to the square lattice at left 
(right). Superconducting reservoirs (with zero superconduct- 
ing phase) are connected on top and bottom. The normal and 
superconducting electrodes are modeled by a collection of one 
dimensional channels. 



CAR transmission in the non-local conductance. We 
also show that an applied magnetic field does not mod- 
ify the balance between EC and CAR. The value of the 
EC and CAR transmission coefficients increases upon ap- 
plication of voltage or magnetic field (positive magneto- 
transmission). For simplicity we do not take into account 
in the present work the effect of Coulomb interactions, 
although it has been realized that they can play an impor- 
tant role for the case of low barrier transparencie s^^' 
We also restrict our analysis to the ballistic case at zero 
temperature. 

In addition to providing insight into the gap and cur- 
rent profiles in the superconducting region, our work 
demonstrates that a description of non-local transport 
in terms of linear relations involving the EC and CAR 
transmission coefficients only makes sense in a quasi- 
equilibrium situation when either the transparency of the 
S/N interfaces or the bias voltage is low. If the supercon- 
ductor is driven far from the equilibrium situation this 
description breaks down, and the resulting non-local con- 
ductance and resistance deviate considerably from their 
values in equilibrium. This case corresponds to the ex- 
perimental situation of Ref. where the non-local resis- 
tance was measured for currents close to the critical one. 
Within our model we are able to obtain the observed 
change of sign of the non-local resistance as a function of 
the injected current^. In particular we demonstrate that 
this change of sign is related to the suppression of the 
local conductance and not due to the dominance of CAR 
over EC processes. 

The article is organized as follows. In the next Sec. |TT] 
we describe the model for a multiterminal NSN structure 
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based on a tight-binding Hamiltonian with a local pair- 
ing which will be determined self-consistently. We also 
give details on how the electronic and transport prop- 
erties of this model are obtained with the help of non- 
equilibrium Green functions. In Sec. Ill we present re- 
sults for the behavior of the complex order parameter and 
the current profiles in the presence of an applied mag- 
netic field, applied voltage and different values of the 
barrier transparency. We also analyze the behavior of 
the non-local conductance for both the quasi-equilibrium 
and non-equilibrium regimes. Concluding remarks are 
provided in Sec. IIVI 



II. MODEL FOR A MULTITERMINAL NSN 
STRUCTURE 



We consider a superconducting region whose thickness 
is much smaller than the coherence and the London pene- 
tration lengths. This region is connected to four electron 
reservoirs as shown in Fig. 2. 

We describe the central superconducting region by 
means of a tight-binding model on a square lattice: 



(1) 



The variable Ak in Eq. ([T]) is the superconducting or- 
der parameter at site k and a summation over pairs of 
neighboring sites {k, I) is carried out in the kinetic term. 
In order to take into account the effect of a magnetic field 
B = rotA one should introduce a phase in the hopping 
elements 



tkl 



texp / A(r).dr 



(2) 



where A(r) is the vector potential, and (j)o = h/e is the 
flux quantum. In the following the magnetic field will be 
measured in terms of (j)/4>o, i.e. the total flux through 
the central region in units of 0o. Disorder could be in- 
troduced in the form of an on-site random potential on 
each tight-binding site. We use the notations Rx = Nao 
and Ry = Mag for the dimensions of the lattice, where 
ao denotes the lattice spacing. 

The central superconducting region is connected at the 
left to the normal electrode and at the right to the 
normal electrode Nq. In order to model a real experi- 
mental situation for non-local transport we connect the 
top and bottom of the central region to superconducting 
reservoirs (see Fig. [2|) by highly transparent interfaces. 
Thus, a current injected through the left interface can 
flow into these reservoirs. The leads, both the normal 
and superconducting, are described by independent one 
dimensional wires connected to each of the lateral sites of 
the central superconducting region. By varying the hop- 
ping terms ta and <{, connecting the central region to the 



normal leads one can control the interface transparency 
Ttv between and 1. The hopping parameter in the rest 
of the system is assumed to have the same value t. 



Calculational Methods 

The transport properties of this model can be conve- 
niently expressed in terms of the retarded G"", advanced 
and Keldysh G"' Green functions, which are matri- 
ces in the Nambu space and depend on two site labels 
k,l. 

The Nambu representation of the central region Hamil- 
tonian takes the form 
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where tki is restricted to first neighbors and was defined 
inEq. 1^. 

The advanced and retarded Green functions are ob- 
tained from a recursive algorithm. To this end we divide 
the central square lattice into layers along the y direction 
which we label by an index n (1 < n < N). We denote by 
G° „ the Nambu advanced Green function projected on 
the sites corresponding to the n layer. This is obtained 
as 



(4) 



— T^n,n+lSn+l'^ri+l,n 

where the recursions take the form 



(5) 
(6) 



with the boundary conditions gg = gjv+i = {i/t)i, 
corresponding to the uncoupled normal wires in the 
wide band approximation. In the recursive equations 
(S„)fc; = hki +taz gs(^zSkiiSki +SkM) with gs = {-luI + 
Aq&x)/ \/ A§ — w^, is the local self-energy on the n layer 
and T„.„±i contains the hopping elements connecting 
neighboring layers. In the last expression we have in- 
troduced the usual Pauli matrices in Nambu space Oz 
and Ox- The dependence on energy u is implicit in the 
Green functions in Eqs. ([I])-®. The Green function con- 
necting two arbitrary layers n and m can then be ob- 
tained from the relations G5^^„j = G° „j_iTm_i^mg^ or 

^n.rn = ,„+iT,„+i,™g^. In the absence of voltage 
applied on electrode N;,, the Keldysh Green function of 
the superconductor takes the equilibrium form 



(7) 



where n(a;) is the Fermi distribution function with zero 
chemical potential. When a voltage Vf, is applied on elec- 
trode Nh the Keldysh Green function of row n is given 
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FIG. 3: (Color online). Self-consistent gap profile at zero volt- 
age for a system size Rx/ao = Ry/ do = 20 and two different 
values of the normal transparency. 



by 



Gt-iuj) = G+-'^«(^) (8) 
+ G«i(c.)Ti,o<Sg;^-(c.)to,iG^„(^)' 



where 



2i ( Sn, 



<5go-H = y, ^^^^^ 



with (assuming zero temperature) 

dueicu) = e{-LO + eVb) - e{-uj) 



(9) 



(10) 
(11) 
I 



The pairing ampUtude in the superconductor is set by 
the anomalous component of the Keldysh Green function 
given by Eq. ([5]). Self-consistency equations at each site 
k take the form 



2iTr 



12 



(12) 



where A is the strength of the attractive electron-electron 
interaction. This is chosen in order to obtain the same 
gap parameter Aq as in the upper and lower supercon- 
ducting leads in equilibrium conditions. Eq. (|12p is iter- 
ated until self-consistency is achieved. The self-consistent 
gap, the current flow in the superconductor and the non- 
local conductance are then evaluated. The current flow- 
ing from site k to site I is given by 



Iki — -r 



duj 



iikGt' 



11 



(13) 



A stringent test of self-consistency is provided by cur- 
rent conservation at each site of the central region be- 
cause current is conserved only once self-consistency is 
verifled^ . 

The non-local conductance can be computed as 

gaAVa,Vb) = dIa/dVb{Va,Vt), whcrC Ia{Va,Vb) = 

J2p-^apap{Va, Vb) is thc total current flowing to electrode 
Na in response to voltages Va and Vb on electrodes a and 
b. We use here Va = 0, as in available expcriments^i^'^. 
The notation Op {bq) is used for the site at thc right (left) 
interface formed by the superconductor and thc normal 
channel p (q), while ap {(5q) are used for the counterpart 
in the superconductor (see Fig. [2]). In quasi-equilibrium 
conditions, i.e. when the variation of the self-consistent 
order parameter with Vb is negligible, the non-local con- 
ductance can be written as^ 



Qa,b{eVb) 



{TcARieVb) - TEcieVb)) 



(14) 



where the EC and CAR transmission coefficients in the 
present model are given by 



TEcieVb) = 2tltl/t'J2{[G^c.p,pMVb) 
TcARieVb) - 2tltl/fY.{[Gk^pS^^^) 



Gl,^^{eVb)\^^ + [Glp,p^eVb)\^^ [Gl^CpieVb) 
Gl^c^pieVb)]^^ + [G^^^^p^ieVb)]^^ [Gl^a^eVb) 



(15) 



(16) 



The sum overp and q in Eqs. (|15p and p6|) corresponds electrodes and N;, to the superconductor, 
to a sum over all one dimensional transverse channels in 
the normal electrodes. As mentioned above, the param- 
eters ta and tb are the hopping amplitudes connecting 
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(a) Tm=0.04, a Vj,/Ao=0.9 (b) 7^=0.84, a Vj,/Ao=0.9 (c) Tfj=0.64, e V[,/A(j=0.5 




(d) Tn=0.04, a Vt,/Ao=0.9 (a) 7^=0.84, a Vj,/Ao=0.9 (f) 1^=0.64, a V|,/Ao=0.5 
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FIG. 4: (Color online). Order parameter and current profiles for a sample size Rx/ao = 12, Ry/ao — 20 under an applied 
voltage Vb- The upper panels show the profiles of the self-consistent gap. The gap is normalized to its valued Ao in the 
superconducting reservoirs. The voltage is indicated by an arrow on the vertical axis. The middle panels show the phase profile 
and the lower ones the corresponding current maps. The size of the arrows in the three lower panels is proportional to the value 
of the local current, with the same scaling factor for the three panels. The values of the parameters Tjv and Vb are indicated 
on top of each panel. 



III. RESULTS 
A. Self-consistent gap, phase and current profiles 

As a first step it is instructive to analyze the behavior 
of the complex order parameter which is obtained from 
the self-consistency Eq. (fT2|) . In order to be able to ac- 
curately describe the spatial variations with a reasonable 
computational cost we fix Ao/t = 0.1, which roughly 
corresponds to a coherence length — 7ao. The system 
sizes that we consider in this work correspond typically 
to Rx,Ry ~ 12 — 30ao but in some cases we use Rx up 
to 60ao and Ry up to 40ao. 

The profile of the self-consistent order parameter am- 
plitude (gap) at zero voltage and zero magnetic field 
is shown on Fig. [3] for a low {Tm — 0.04) and a high 
(Tat = 0.64) values of the interface transparencies. As 
can be observed the gap is reduced at the contacts with 
the normal leads due to the inverse proximity effect. This 
gap suppression becomes stronger as the interface trans- 
parency increases. There is also a slight suppression on 
the contacts with the upper and lower superconducting 
reservoirs which appears due to its deviation from an 



ideal interface. On the other hand, in equilibrium and 
for sufficiently large system size the gap in the middle of 
the central superconducting region is close to its value in 
the reservoirs. One can also notice that the gap profiles 
exhibit small ripples along certain directions, which in 
the case of a square region {R^ = Ry) correspond to the 
diagonal lines. These ripples reflect the ballistic charac- 
ter of our model which leads to constructive interference 
along semi-classical trajectories (for a discussion see Ap- 
pendix) . 

An applied voltage on the left electrode N;, leads to sev- 
eral modifications in the order parameter profiles. This 
is illustrated in Figs. |4] and [5] for system sizes R^ = 12ao 
and Rx = 30ao- First, there is a reduction in the aver- 
age gap value due to dcpairing effects. This reduction 
is more pronounced near the contact where the current 
is injected. It is interesting to note that for the short 
system and high transparency [panel (b) in Fig. [4] the 
average gap becomes smaller than the applied voltage for 
Vf,/Ao = 0.9. In contrast for samples with Rx sufficiently 
large compared to the coherence length, the overall gap 
profile is less sensitive to the applied voltage even for high 
transparency [panel (b) in Fig. [S]. 
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(b) T,=0.64. e V,/A„=0.5 1^=0.64, 6 Vt,/Ao=0.5 






(b) Tm=0.64, e \l&o=0.5 



(e) 1^.0.04, e V|,/A„=0.5 



(f) T„=0.64, e V|,/A„=0.5 
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FIG. 5: (Color online). Same as Fig. |4]for sample dimensions 
Rx/ao = 30 and Ry/ao — 20. The scaling factor between 
current and size of arrow is not the same as in Figs. |4^, h and 
i, and it is different for panels e and f. 



Secondly, the application of a finite voltage Vb leads to 
a smooth drop of the phase of the order parameter along 
the X direction [panels (d)-(f) in Fig. [4] and (c)-(d) in 
Fig. As can be observed this is accompanied by a 
phase gradient in the y direction which changes sign at 
the midpoint y ~ Ry/'2. 

The properties of the self-consistent solution can be 
understood more physically by analyzing the current pro- 
files. These arc shown in the lower panels of Figs. |4]and 
[51 As a general remark, the profiles illustrate how the in- 
jected current from the left electrode gradually leaks to 
the upper and lower superconducting electrodes. As ex- 
pected, the amount of current reaching the right normal 
lead is controlled by the barrier transparency, the applied 
voltage and the system size. For sufficiently large system 
size one can clearly appreciate an exponential suppres- 
sion of the injected current on the scale [panels (e) 
and (f) in Fig. In addition the current profiles in the 
case of the larger system exhibit some structure along 
the lines y ~ x and y = Ry — x which can be associated 
with the already mentioned ripples appearing in the gap 
profiles. 

In the presence of applied voltage and magnetic field 
the phase profile becomes more complex with a modula- 
tion both in the x and y directions. As shown in Fig. [6] 
this is associated with the appearance of currents along 
the y direction which tend to screen the applied field. 
These currents are essentially superimposed to the pro- 
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FIG. 6; (Color online), (a) Surface map of the self-consistent 
phase under magnetic field, for sample dimensions Rx/ao — 
12 and Ry/ao = 20. The corresponding current map is shown 
on (b). The flux is (j}/(j)o = 1. The scaling factor between 
current and size of arrow is the same as in Figs. [4^, h and i. 



file arising from the injected current through the left elec- 
trode. On the other hand, for this range of parameters 
the effect of the magnetic field on the amplitude of the 
order parameter corresponds to a uniform reduction (not 
shown here). 

The superposition of injected and screening currents 
can be clearly observed in the evolution of the current 
profiles for increasing values of the magnetic field, shown 
in Figs. [7] and [8] for the cases Rx = 12ao and Rx = 
30ao respectively. These plots correspond to the case of 
high transparency (T/v = 0.64) and high voltage {Vi, = 
0.9Ao) for which the injected and screening currents have 
a similar size when (f'/'Po — 1- For a given total flux (j) 
the screening currents are, however, smaller for the case 
of the shorter system due to the larger gap suppression 
induced by size effects and by the applied voltage. 



B. Low bias non-local conductance 

In this subsection we analyze the non-local transport 
properties in the regime eVj, <^ Aq. As discussed in the 
section on calculational methods the non-local conduc- 
tance in quasi-equilibrium conditions can be decomposed 
into EC and CAR contributions, given by Eqs. (fT5)l and 
(fT6|l . These two contributions cannot be disentangled 
in experiments, but theoretically it is convenient to ana- 
lyze them in a separate way, especially in the case of low 
transparency tunnel barriers where the two contributions 
nearly cancel in the non-local conductancei. 

The variation of the low bias EC and CAR transmis- 
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FIG. 7: (Color online). Current map for increasing values of the magnetic flux 0/0o and for T/v = 0.64, Rx/ao = 12 and 
Ry/ao = 20. The applied voltage 14 is set to 0.9Ao. An arrow with the same length corresponds to the same current for all 
panels. The scaling factor between current and size of arrow is the same as in Figs. [4^, h and i. 
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FIG. 8: (Color online). Same as in Fig.[7]for Tjv = 0.64, eVb = 0.9Ao, R^/ao = 30 and Ry/ao = 20. 



sion coefficients with barrier transparency Tn is shown 
on Fig.ini We find it convenient to represent in this figure 
the quantities Fec^Tn) and Fcar{Ti^) defined as 



Fec{Tn) 
Fcar{Tn) 



{Tn)-^Tec{Tn) 
hmT„ -^o{TNy^TEc{TN) 

{Tn)-^Tcar(Tn) 
lh-nTf,^oiTN)~'^TEc{TN) ' 



(17) 
(18) 



This normahzation allows to analyze the variation of 
Tec a-nd Tqar beyond the dominant ~ Tfj dependence 
which is common to both coefficients. As can be ob- 
served, transmission of electrons dominates over crossed 
Andreev processes for highly transparent interfaces while 
they become almost equal in the T/v — > limit. This is 
an effect previously discussed in literature^i^ which is not 
modified significantly in our self-consistent calculation. 

It is also interesting to analyze the dependence of the 
EC and CAR coefficients on the superconductor length 



Rx- This is illustrated in Fig.fTOl We choose to represent 
in this figure only Tec as Tcar. exhibits the same dis- 
tance behavior. We find that except for tiny fluctuations 
at certain values of Rx the normalized transmission coef- 
ficients follow an exponential decay with sample dimen- 
sions of the type T{Rx) ~ exp (— 2i?a;/^e//), where ^e// 
is an effective coherence length which is sensitive to the 
applied voltage and magnetic field. As shown in panel (a) 
of Fig. [TOl the effective coherence length is also sensitive 
to the transverse dimension Ry , increasing slightly as Ry 
increases. From Fig. [TOk we can estimate the values of 
^e//. The coherence length for Ry = 15ao is ^e// ~ 7ao, 
while for Ry = 25ao is S.eff ~ S.Soq. The increase of ^e// 
by increasing Ry is due to the enhancement of the inverse 
proximity effect which reduces the average gap in the S 
region. 

On the other hand, panel (b) on Fig. [TUl illustrates the 
increase of ^eff when a magnetic field is applied. The 
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FIG. 9: (Color online). Variation of the quantities Fec(Tn) 
and Fcar{Tn) defined in text, as a function of the normal 
transparency Tjv. A predominance of EC over CAR is ob- 
tained for high values of interface transparencies. The chosen 
system dimensions are Rx = 28ao and Ry — 20ao. 
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FIG. 11: (Color online). Dependence of the transmission co- 
efficients Tec {(I}) /Tec {(t> = 0), and TcAR{4')/TEc{<j} = 0) 
on the magnetic flux (j>/4'o- We have chosen Tn = 0.64, 
Rx/a.0 = 20 and Ry/ao = 14. The broken and dotted lines 
indicate the quadratic fit of Ea. p9p . 



Finally, Fig. [11] shows the flux dependence of the EC 
and CAR transmission coefficients. We find that this 
dependence can be fitted as 



Tec,car{4>) = Tec,car{4' = 0) 



0* 



,(19) 



6 8 10 

2Rx/^eff(0) 



where 4)*/4>{) is roughly of order unity. For instance, the 
same value (t)*/<Po ~ 1.35 ± 0.05 for the EC and CAR 
transmission coefficient is obtained within error-bars for 
the data on Fig. [TTJ 



C. Non-local conductance at arbitrary bias 



FIG. 10: (Color online). Dependence of the EC transmission 
coefficient on the length Rx = A'^ao of the superconductor. 
Panel (a) corresponds to low transparency and two different 
values of the transverse dimension Ry = 15 and 25ao. The 
sample dimension Rx on panel (a) is normalized to the value 
Ce//(0) of the coherence length for Ry — 15ao. Panel (b) 
corresponds to high transparency and two diff'erent values of 
the magnetic flux. The sample dimension on panel (b) is 
normalized to the value Ce//(0) of the coherence length for 
(f>/4>0 = 0. Out-of-equilibrium effects are negligible for both 
panels (see insets of Fig. I12p . 

coherence length for 0/00 — 0.6 is 20 % larger than for 
4>/(t>o = 0, which again can be associated to the reduction 
of the average gap, now due to the depairing effect of the 
applied magnetic field. 

As can be more clearly observed for low transparency, 
the transmission coefficients show small fluctuations 
around ~ pRy, with p being an integer. This is an- 
other feature associated with the ballistic character of 
our model as it is discussed in the Appendix. 



At arbitrarily large bias the non-local conductance can- 
not be obtained from Eq. (jl4p with the transmission co- 
efficients defined by Eqs. and p6|) as it contains 
contributions due to the voltage dependence of the self- 
consistent order parameter in the S region. More gener- 
ally it can be computed directly as dIa/dVb- 

The upper panels of Fig. [12] illustrate the behavior 
of Qa,b at arbitrary bias for a system size {Rx,Ry) = 
(14ao, 20ao) and two different values of the normal trans- 
parency Tn. For small Tn the non-local conductance ex- 
hibits a abrupt jump for eVf, Aq. In this limit Qa^b is 
well described by Eq. (fT4|) [see inset in panel (a)] which 
is a consequence of having a quasi-equilibrium situation 
in the whole voltage range. This is also reflected in the 
magnetic field dependence which remains similar to the 
one found at low bias. In contrast, for high transparency 
non-equilibrium effects manifest in several features of the 
non-local conductance. For instance one can clearly no- 
tice that the jump associated to the gap appears at lower 
voltage bias (Fig[T2]b). Also, as shown in the inset, the 
actual value of Qa.b deviates from the one calculated from 



9 



(a) Tn = 4 10- 



— 


- 










nj 










-5 - 




















-10 - 






Clf 




(M 




o 






-15 - 









100 r 








75 - 








50 - 






> 




3- 


25 - 












- 




(b) Tn = 0.64 






■i - 








2.5 - 


OC 




< 


2 - 


> 

3, 


1.5 - 








1 - 




0.5 1 1.5 

eV^/Ao 



FIG. 12: (Color online). Voltage dependence of the non-local conductance (upper panels), and of non-local resistance (lower 
panels), normalized to their value at zero bias in zero applied magnetic field. The values of Tn and (j) are indicated on each 
panel. The insets of the upper panels illustrate the comparison between the full calculation Qa,b ~ dla/dVt (solid line) and the 
expression approximation Qa,b = /h{TcAR — Tec) with the transmission coefficients defined by Eqs. p5[) and p6|) (dashed 
line). The non local conductance Qa,b in the insets is given in units of e^//i. 



(a) Rx=13ao 



(b) R,=15ao 



a -250 
-500 

( 

0.1 

„ 0.05 

£ 

-0.05 
-0.1 




0.2 0.3 0.4 0.5 0.6 0.7 

e Vi,/A(, 

(c) R,=13ao 




0.2 0.3 0.4 0.5 0.6 0.7 

e \IAo 



a -250 
-500 

( 

0.1 

_ 0.05 

CD 

£ 

cc" 

-0.05 
-0.1 




0.2 0.3 0.4 0.5 0.6 

e V^/Aq 

(d) R,=15ao 







Ra,b , 1 





0.2 0.3 0.4 0.5 0.6 

e Vi,/At, 



0.7 O.E 



FIG. 13; (Color online). Voltage dependence of the local and non-local conductances (upper panels) and of the non-local 
resistance (lower panels) for an asymmetric junction with Tjv = 1 at the b interface and Tjv = 0.1 at the a side, for Ry = 40ao 
and two values of R^: 13ao [panels (a) and (b)] and 15ao [panels (c) and (d)]. 



Eq. ([T4ll . Finally, it is found that the magnetic field at 
voltages ^ Aq reduces the non-local conductance, i.e. it 
has the opposite effect to the one found at low bias. 

In order to make contact with existing experiments it 
is also interesting to analyze the behavior of the non-local 



resistance, defined as 

Tla.b = -Ga.b/[Qa,aGb.b — Ga,bGb,a] (20) 

where Ga,a and Gb,b denote the local conductances at 
each interface and Ga,b and Gb,a are the non local con- 
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ductances. The later take approximately the same value 
under quasi-equilibrium conditions, but they are in gen- 
eral different due to the voltage dependence of the super- 
conducting order parameter. Eq. (j20p is valid provided 
that the induced voltage Va is very small and under the 
assumption that the current voltage relations can be lin- 
earized around the voltage Vb under consideration. The 
corresponding results are shown in the lower panels of 
Fig. [121 While TZa.b is almost constant in the region 
of small voltages it exhibits a high increase around the 
self-consistent gap. 

Even stronger non-equilibrium effects arc found when 
increasing the transverse dimension Ry and for higher 
values of the normal transparency Tn at the b side while 
maintaining at a at low values. In that case most of the 
injected current flows into the superconducting electrode 
and one can obtain an abrupt jump in TZa,b from posi- 
tive to negative for Vb of the order of the self-consistent 
gap>S. The origin of this change of sign can be under- 
stood from the results in Fig. [13] where the local and 
non-local conductances are plotted for different values of 
Rx (upper panels). As can be observed, the local con- 
ductances exhibit a decrease for voltages approaching the 
self-consistent gap. This decrease can even lead to a neg- 
ative local differential conductance which arises due to 
the fact that the injected current can become larger than 
the maximum supercurrent that can leak into the super- 
conducting electrode. A similar effect was described in 
Ref. [33. This effect is reduced by increasing the length 

As a consequence of the suppression of the local con- 
ductances there exists a voltage window where one can 
have ga,bQb,a > Ga,aQb,b which, according to Eq ([2Ql), 
leads to a change of sign in the non-local resistance. It is 
worth emphasizing that this change of sign is not related 
to a change of sign of Qa^b and therefore not associated to 
a possible dominance of CAR over EC processes. These 
results arc in qualitative agreement with the experimen- 
tal data of Ref. [^. 

We plot in the lower panels of Fig. [TSjthe variation of 
TZa,b as a function of Vb- One can observe that the abrupt 
change of sign in TZa,b reduces its amplitude and shifts 
towards higher voltages as Rx increases. The shift corre- 
sponds to the increase of the critical current of the central 
region which grows linearly with R^, an effect which is 
not present in the data of Ref. [^ as it corresponds to a 
somewhat different geometry. Moreover, while our calcu- 
lations have been performed in the voltage biased case, 
the experiments of Ref. correspond to a current biased 
situation. Thus, the regions of negative differential con- 
ductance could not be explored. However, the regions 
with negative TZa,b shown in the lower panels of Fig. [T3| 
would be observable, since the sign change of TZa,b is re- 
lated to the condition Ga.aGb.b < Qa,bQb.a rather than to 
the change of sign of Qb,b. In fact, we have checked that 
the change of sign of TZa.b takes place outside the regions 
with negative Gb,b- 
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FIG. 14: (Color online). Dependence of the normalized EC 
and CAR transmission coefficients T{Rx)/T{0) on 2Rx/^q, 
where Rx is the dimension of the superconductor along x axis, 
and 5o is the bulk coherence length. The applied voltage is 
small compared to the gap {eVt/Ao — 0.02). The coherence 
length is larger than for the corresponding parameters that in 
the presence of superconducting leads. Panel (a) corresponds 
to Ry — 5ao and panel (b) to Ry — 25 ao. The symbols 
correspond to EC and the bold red lines to CAR. 



IV. CONCLUSIONS 

We have presented self-consistent model calculations 
for describing a typical experimental setup to measure 
non-local transport in nanoscale superconductors con- 
nected to normal electrodes. An important ingredient in 
these calculations is the inclusion of the superconducting 
leads which drain the current injected into the nanoscale 
region. 

We have found that this system exhibits a rich vari- 
ety of behaviors controlled by several parameters. One 
can distinguish two limiting cases: the quasi-equilibrium 
regime in which the superconducting order parameter 
depends only weakly on the applied voltage and the 
strong non-equilibrium case where this dependence can- 
not be neglected. In these two regimes the non-local 
conductance exhibits a quite different behavior. In quasi- 
equilibrium conditions Qa,b can be expressed as the differ- 
ence between Tcar and Tec- We have shown that Tec 
becomes larger than Tcar for increasing transparency, 
and that both decay exponentially with the system size 
Rx- Moreover we have shown that both coefficients in- 
crease under an applied magnetic field. 

In the case of strong non-equilibrium a simple expres- 
sion of Ga,b as the difference between Tcar and Tec is 
no longer possible. The voltage dependence of the self- 
consistent order parameter introduces an additional con- 
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tribution to the non-local conductance. We have also 
found that even the local conductance can be strongly 
modified by non-equilibrium eS'ects. The efi^ects on these 
two quantities lead to a non-local resistance which may 
exhibit a strongly non-monotonous behavior as a func- 
tion of the injected current including a sign change for 
sufficiently large transparency. 

As a final remark we would like to comment the pos- 
sible connection between the present model calculations 
and the existing experiments on non-local transport. As 
mentioned in the introduction, basically three different 
experiments have been performed so far. They corre- 
spond to different geometries and to different materi- 
als. The size of the junctions transparency can also 
change enormously from one experiment to the other. 
It is thus not possible to infer general conclusions from 
these results. In the case of the experiment of Rcf. 
the authors claim to have junctions with a extremely 
small transparency ^ 10^^ which would warrant the 
quasi-equilibrium conditions in the subgap voltage range. 
The results of the present work, corresponding to a non- 
interacting theory, would not be able to describe that 
experiment. In particular our results predict an increase 
of the non-local conductance with magnetic field in con- 
tradiction with the experimental observations. As stated 
in a previous work by some of the authors^, interactions 
mediated by the electromagnetic modes may play an im- 
portant role in that experimental situation. It was also 
conjecturccU^ that disorder may play a role when associ- 
ated to phase fluctuations. 

On the other hand, the experimental results of Ref. 
clearly correspond to the strong non-equilibrium situa- 
tion described in the present work. Even when our model 
geometry does not fully correspond the their experimen- 
tal setup the behavior of the non-local resistance as a 
function of the injected current shown in Fig [13] is very 
similar to the one found in this experiment. A closer 

I 



[5a,/3(w)]n 



where the quasiparticle energy En^m is such that E^^ = 

A' + iUm?, with Cn.™ = 2f[c0s(^) + C0s(^)]. 

The coherence factors are given by {un,m)^ = (1/2) (1 + 

^n,7n / Eji_jji^ ^ (^n.m) — (1-/2)(1 Cn,m/-^n,m)- The 

anomalous component [(7q./3]j^2 obtained by substitut- 
ing the numerator uj+iri+^n,m in Eq. ((AT|) by - A. These 
Green functions lead to an enhanced probability for the 
propagation along certain directions. As one would ex- 
pect from a semi-classical analysis of an integrable cav- 



comparison between theory and experiments would be 
desirable for further understanding of the observed fea- 
tures. 
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APPENDIX: GREEN FUNCTIONS AND 
TRANSMISSION COEFFICIENTS OF A CLOSED 
RECTANGULAR SUPERCONDUCTING CAVITY 

In order to understand the oscillations in the EC and 
CAR transmission coefficients as a function of the length 
Rx (Fig llOp . we consider in this Appendix a simpli- 
fied model consisting of a central superconducting re- 
gion which we describe by the tight-binding Hamilto- 
nian of Eq. ([T|) on a sample of dimensions x Ry, with 
Rx — Nqq and Ry = Muq. We do not include the con- 
tacts to the superconducting leads and do not implement 
self-consistency for the order parameter. 

We assume tunnel contacts with the normal electrodes 
at left and right so that the non-local conductance is 
given by lowest order perturbation theory in the tunnel 
amplitudes. The normal Green function describing a par- 
ticle propagating from site a at coordinates {xa,ya) to 
site (3 at coordinates {xp,yp) is given 



(A.l) 

I 

ity these directions correspond to the diagonal lines in 
the case of a square superconducting region. The rip- 
ples shown by the self-consistent gap are a consequence 
of constructive interference along these semi-classical tra- 
jectories. 

The EC and CAR transmission coefficients can be eval- 
uated using these approximated Green functions. Fig. 1141 
shows the behavior of the transmission coefficients T^c 
and TcAR as a function of Rx obtained with this sim- 



= (^) (^) V V sin + M sin f„ 

Vl + Af/Vl + ^^/n, V N+1 J \ N+l 

. ( ya/ao + l\ . / 2/^3/00 -I- l\r t^ + M7 + ?ri.m 

sm m — — ; — sm m^-— ; — 7 — -rr. — -7- tt: 



12 



plified model for two values of Ry. We find that the 
the EC and CAR transmission coefficients almost coin- 
cide for Rx in between pRy and {p + l)Ry, where p is an 
integer. Pronounced oscillations are obtained when the 
dimensions of the superconducting region are such that 
Rx/Ry — P- Again these effects are a consequence of 



interferences along semi-classical trajectories in this bal- 
listic intcgrable model. In the same way, one could use 
these approximate Green functions in a self-consistent 
calculation to show that the ripples in the self-consistent 
gap (see e.g. Fig.^ have a similar origin. 
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